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Abstract 

Detailed models of Procyon A based on new asteroseismic measurements by Eggen- 
berger et al. (2004b) have been computed using the Geneva evolution code including 
shellular rotation and atomic diffusion. By combining all non-asteroseismic observ- 
ables now available for Procyon A with these seismological data, we find that the 
CN ■ observed mean large spacing of 55.5 ± 0.5 /^Hz favours a mass of 1.497 M@ for Pro- 

cyon A. We also determine the following global parameters of Procyon A: an age 
\ of t = 1.72 ± 0.30 Gyr, an initial helium mass fraction Yi = 0.290 db 0.010, a nearly 

solar initial metallicity (Z/X)\ = 0.0234 ± 0.0015 and a mixing-length parameter 
a = 1.75 ± 0.40. Moreover, we show that the effects of rotation on the inner struc- 
ture of the star may be revealed by asteroseismic observations if frequencies can be 
determined with a high precision. Existing seismological data of Procyon A are un- 
fortunately not accurate enough to really test these differences in the input physics 
of our models. 
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1 Introduction 



Due to its brightness and proximity, Procyon A constitutes an ideal target to 
test the input physics of the stellar models and to search for p-mode oscilla- 
tions. 

Hartmann et al. (1975) were the first to calculate stellar models of Procyon A 
and to compare them to the observational data available at the time. They 
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found a discrepancy between the astrometric mass (Strand 1951) and the 
lower astrophysical mass deduced from their models. This mass discrepancy 
was later confirmed by Demarque & Guenther (1988), who also found that 
overshoot of the convective core had to be included in the models in order to 
reproduce the observed effective temperature. 

To adress this problem of mass discrepancy, Irwin et al. (1992) decided to 
redetermine the astrometric mass of Procyon A. They found a mass in perfect 
agreement with the previous value of 1.74 M Q . Following this redetermination 
of Procyon's mass and motivated by several attempts made to detect the 
signature of oscillation modes, Guenther & Demarque (1993) calculated an 
array of stellar models for Procyon A including the new OPAL opacities. They 
concluded that, using OPAL opacities, no convective overshoot was needed 
to match Procyon's position in the HR diagram. However, the discrepancy 
between the astrometric mass and the mass deduced from stellar evolution 
models still remained. 

This discrepancy was solved thanks to new measurements of orbital elements 
and parallax of the Procyon system by Girard et al. (1996, 2000). They de- 
termined an astrometric mass of 1.5 M Q in good agreement with the value 
supported by stellar models. Using this revised astrometric mass, Chaboyer 
et al. (1999) calculated a grid of stellar evolution models for Procyon A and 
investigated their seismic properties. They concluded that the detection of 
p-modes would serve as a robust test of stellar evolution theory. 

The first indication of the presence of p-modes on Procyon A was obtained 
by Brown et al. (1991), while the first clear detection was made by Martic et 
al. (1999). By comparing these observations with theoretical predictions and 
numerical simulations, Barban et al. (1999) confirmed the stellar origin of the 
observed excess power. 

Very recently, individual p-mode frequencies were identified (Eggenberger et 
al. 2004b, hereafter ECBB04; Martic et al. 2004). However, Matthews et al. 
(2004) reported that the Most satellite did not observed any evidence of the 
expected acoustic oscillations. To explain the apparent discrepancy with previ- 
ous radial- velocity measurements, Matthews et al. suggested that either both 
the radial velocity measurements and the Most observations are dominated 
by granulation noise, or the properties of the oscillations are different from 
the theoretical expectations. Christensen-Dalsgaard & Kjeldsen (2004) cast 
doubt on the fact that granulation dominates the noise, and suspected that 
the Most data might be dominated by non-stellar noise. Moreover, new as- 
teroseismic observations with the Harps spectrograph confirmed the previous 
Doppler ground-based detections (Bouchy et al. 2004). 

In this work, we will combine all non-asteroseismic measurements with aster- 
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oseismic observations of ECBB04 to investigate which additional constraints 
are brought by these seismological data. We will thus try to determine a model 
of Procyon A which best reproduces all these observational constraints using 
the Geneva evolution code which includes a complete treatment of shellular 
rotation and atomic diffusion. Moreover, we will investigate the effects of ro- 
tation on the global parameters of Procyon A and on the p-mode frequencies. 

The observational constraints available for Procyon A are summarized in 
Sect. 2, while the input physics of the models and the calibration method 
are described in Sect. 3. The results are presented in Sect. 4 and the conclu- 
sion is given in Sect. 5. 



2 Observational constraints 

2.1 Astrometric data 



The astrometric parameters of the visual binary orbit were recently updated by 
Girard et al. (2000). Using data obtained with the infrared cold coronagraph 
(CoCo) and with the WFPC2 on board the Hubble Space Telescope, they 
found a mass of 1.497 ± 0.037 M & for Procyon A, with a parallax IT = 283.2 ± 
1.5 mas. When only the WFPC2 observations are used, a smaller mass is found: 
M = 1.465±0.041 M Q . These masses were derived by using the self-consistent 
parallax of 283.2 ± 1.5 mas. Using the slightly larger parallax measured by 
Hipparcos (n = 285.93 ± 0.88 mas) with the astrometric parameters of the 
orbit results in a smaller mass than the one obtained with the parallax of 
Girard et al. (2000). Indeed, when only the WFPC2 measurements are used, 
we find a mass of 1.423 ± 0.040 M© with the Hipparcos parallax. Thus, we 
see that the mass of Procyon A lies between 1.42 and 1.5 M Q . All of these 
masses will be considered in our analysis. We will then determine which value 
is in best accordance with the other observational constraints and in particular 
with the asteroseismic measurements. 



2.2 Effective temperature and chemical composition 

For the effective temperature of Procyon A, we adopted T c g = 6530 ± 90 K 
(Fuhrmann et al. 1997). Note that the value of 6530 K is also given by Allende 
Prieto et al. (2002). The chemical composition of Procyon A is nearly solar; 
we adopted [Fe/H] = -0.05 ± 0.03 (Allende Prieto et al. 2002). 
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2.3 Luminosity 



From the compilation of 13 measurements from the literature, Allende Prieto 
et al. (2002) derived a mean visual magnitude (V) = 0.363 ± 0.003 mag for 
Procyon A. Combining this mean magnitude with the solar absolute bolo- 
metric magnitude Mb i,© = 4.746 (Lejeune et al. 1998) and the bolometric 
correction from Flower (1996), we find a luminosity logL/L = 0.84 ± 0.02. 
Note that this interval in luminosity is compatible with the use of the parallax 
of Girard et al. (2000) as well as the Hipparcos parallax. Indeed, the changes 
in L due to the different parallaxes are very small and entirely included in the 
error bar of 0.02 (logL/L Q = 0.843 and 0.835 with the parallax of Girard et 
al. and the Hipparcos parallax respectively). 

2.4 Angular diameter 

Recently Kervella et al. (2004) measured the angular diameter of Procyon A 
using the VINCI instrument installed at ESO's VLT Interferometer. They 
found a limb darkened angular diameter 9 = 5.448 ± 0.053 mas. Using the 
Hipparcos parallax, they deduced a linear diameter of 2.048 ± 0.025 D Q . The 
parallax of Girard et al. (2000) gives a linear diameter of 2.067 ± 0.028 D Q . 

2.5 Rotational velocity 

Allende Prieto et al. (2002) estimated v sini = 3.16 ± 0.50 km s -1 . However, 
they pointed out that the correct value is probably close to 2.7 km s -1 , since the 
value of 3.16 may be slightly overestimated as a result of the finite numerical 
resolution of their convection simulation. Thus, we adopted v sini = 2.7 km s -1 
with a large error of lkms -1 in order to encompass the value of vsini = 
3.16 ± 0.50 kms -1 . Using vsini = 2.7 ± 1.0 kms -1 and assuming that the 
rotation axis of Procyon A is perpendicular to the plane of the visual orbit 
{i = 31.1 ± 0.6°, Girard et al. 2000), we find a surface rotational velocity of 
5.2 ± 1.9 km s -1 for Procyon A. 

2.6 Asteroseismic constraints 

Martic et al. (2004) identified individual frequencies in the power spectrum 
between 300 and 1400 /iHz with a mean large spacing of 53.6 ± 0.5 /zHz. Using 
the Coralie spectrograph, we confirmed the detection of p-modes on Pro- 
cyon A and identified individual frequencies with a slightly larger mean large 
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Table 1 

Observational constraints for Procyon A. References: (1) Girard et al. (2000), (2) 
Hipparcos, (3) Allende Prieto et al. (2002), (4) derived from the other observational 
measurements (see text), (5) Fuhrmann et al. (1997), (6) Kervella et al. (2004), (7) 
Martic et al. (2004) and (8) Eggenberger et al. (2004b). 



References 
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LOO. 1 ± 1.0 
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nor r\o i i \ oo 

285. 9d ± 0.88 


(2) 


M/Mq 


1.497 ±0.037 


(1) 




1.465 ±0.041 


(1) 




1.423 ±0.040 


(l)+(2) 


V [mag] 


0.363 ± 0.003 


(3) 


log L/L Q 


0.84 ± 0.02 


(4) 


T cS [K] 


6530 ± 90 


(5) 


[Fe/H] s 


-0.05 ±0.03 


(3) 


[mas] 


5.448 ± 0.053 


(6) 


R/Rq 


2.048 ± 0.025 


(6)+(2) 




2.067 ± 0.028 


(6)+(l) 


V s [kms- 1 ] 


5.2 ± 1.9 


(4) 




53.6 ±0.5 


(7) 




55.5 ±0.5 


(8) 
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spacing of 55.5±0.5 /iHz (ECBB04). In this work, we will use the asteroseismic 
observations of ECBB04 to constrain the models. 

All the observational constraints are listed in Table 1. 



3 Stellar models 

3.1 Input Physics 

The stellar evolution code used for these computations is the Geneva code in- 
cluding shellular rotation, described several times in the literature (see Meynet 
& Maeder 2000). We used the new horizontal turbulence prescription of Maeder 
(2003) and the braking law of Kawaler (1988) in order to reproduce the mag- 
netic braking that undergo the low mass stars when arriving on the main 
sequence. Two parameters enter this braking law: the saturation velocity fi sat 
and the braking constant K. Following Bouvier et al. (1997), f2 sat was fixed 
to 14 Vl Q and the braking constant K was calibrated on the sun. We used the 
OPAL opacities, the NACRE nuclear reaction rates (Angulo et al. 1999) and 
the standard mixing-length formalism for convection. 

In addition to shellular rotation, our models have been computed including 
atomic diffusion on He, C, N, O, Ne and Mg using the routines developed 
for the Geneva-Toulouse version of our code (see for example Richard et al. 
1996) recently updated by O. Richard (private communication). The diffusion 
coefficients are computed with the prescription by Paquette et al. (1986). We 
included the diffusion due to the concentration and thermal gradients, but the 
radiative acceleration was neglected. 

3.2 Computational method 

The computation of a stellar model for a given star consists in finding the set 
of stellar modeling parameters which best reproduces all observational data 
available for this star. For a given stellar mass, the characteristics of a stellar 
model including the effects of rotation (luminosity, effective temperature, sur- 
face metallicity, surface velocity, frequencies of oscillation modes, etc.) depend 
on five modeling parameters: the age of the star (t hereafter), the mixing- 
length parameter a = l/H p for convection, the initial surface velocity V\ and 
two parameters describing the initial chemical composition of the star. For 
these two parameters, we chose the initial helium abundance Y\ and the ini- 
tial ratio between the mass fraction of heavy elements and hydrogen (Z/X)i. 
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Fig. 1. Observational constraints in the HR diagram for Procyon A. The contin- 
uous lines indicate the boxes delimited by the observed luminosity and effective 
temperature (with their respective 1-sigma errors). The dashed and dotted lines 
denote the boxes delimited by the interferometric radii deduced using the parallax 
of Girard et al. (2000) and the Hipparcos parallax, respectively. The positions in 
the HR diagram of the computed models are indicated (see text and Table 2 for 
more details). 

Assuming that this ratio is proportional to the abundance ratio [Fe/H], we 
can directly relate (Z/X) to [Fe/H] by using the solar value (Z/X) Q = 0.0230 
given by Grevesse & Sauval (1998). Thus, any characteristic A of a given 
stellar model has the following formal dependences with respect to modeling 
parameters : A = A(t, a, V i: Y h (Z/X){). 



Once fixing the mass of Procyon A to one of the three values listed in Table 1, 
the determination of the set of modeling parameters (t, a, V[, YJ, (Z/X){) 
leading to the best agreement with the observational constraints is made in 
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two steps. First, we construct a grid of models with position in the HR dia- 
gram in agreement with the observational values of the luminosity, effective 
temperature and radius listed in Table 1. The boxes in the HR diagram for 
Procyon A corresponding to the three different choices of mass are shown in 
Fig. 1. 

To construct this grid, we proceed in the following way: for a given chemi- 
cal composition (i.e. a given set Yj, (Z/X)i) and a given initial velocity Vi, 
the mixing-length coefficient of each star is adjusted in order to match the 
observational position in the HR diagram. 

Once the position of the star in the HR diagram agrees with the observed 
values, the surface metallicity of the star is compared to the observed one. If 
it is out of the metallicity interval listed in Table 1, the models are rejected 
and the procedure is repeated with another choice of Yj and {ZjX\. Note that 
the surface metallicities [Fe/H] s are almost identical for the models with the 
same initial composition and different mixing-length parameters. Moreover, 
the [Fe/H] s of the models are mainly sensitive to [ZjX)\ and less to Yj. As a 
result, the values of {ZjX\ are directly constrained by the observed surface 
metallicity. 

Finally, the surface velocity of the models is compared to the observed one; if 
the velocities of the models are not compatible with the observed value, the 
models are rejected and the procedure is repeated with another initial velocity 
(but with the same initial chemical composition). Otherwise, all solutions are 
kept, since they correspond to models of Procyon A which reproduce all the 
non-asteroseismic constraints. The whole procedure is then repeated with a 
new choice of Y\ and (Z / X)\. 

In this way we obtain a grid of models with various sets of modeling parameters 
(t, a, Vi, Y i? (Z/X)i) which satisfy all the non-asteroseismic observational 
constraints of Procyon A, namely the effective temperature, the luminosity, 
the radius, the surface velocity and the surface metallicity. The second step 
in determining the best model of Procyon A is to consider the asteroseismic 
measurements. 

For each stellar model of the grid constructed as explained above, low-/ p- 
mode frequencies are calculated using the Aarhus adiabatic pulsation package 
(Christensen-Dalsgaard 1997). Following our observations, modes of degree 
/ < 2 with frequencies between 0.6 and 1.4 mHz are computed and the large 
spacing Az/ determined. 

Once the asteroseismic characteristics of all relevant models are computed, we 
perform a x 2 minimization in order to deduce the set of parameters (t, a, Vi, 
Yi, (Z/X)i) leading to the best agreement with the observations of Procyon 
A. For this purpose, we define the x 2 functional 
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6 / /^theo s^obs \ 2 

where the vectors C contains all the observables for the star: 

C = (log L/L Q , T eff , R/R Q , V s , [Fe/H] s , Au ) . 

The vector C theo contains the theoretical values of these observables for the 
model to be tested, while the values of C obs are those listed in Table 1. The 
vector aC contains the errors on these observations which are also given in 
Table 1. 



4 Results 

4-1 Models with a mass of 1.497 M & 

We first calculated a grid of models including shellular rotation and atomic 
diffusion with the mass of 1.497 M Q deduced from data obtained with the 
infrared cold coronagraph (CoCo) and with the WFPC2 on board the Hubble 
Space Telescope Girard et al. (2000). This mass has been determined with 
a parallax of 283.2 ± 1.5 mas, leading to a radius of 2.067 ± 0.028 R Q for 
Procyon A. The models have to match the location of Procyon A in the HR 
diagram which is given by its radius, luminosity and effective temperature (see 
Fig. 1). 

Once this grid of models was computed, we performed the y 2 minimization 
described above. In this way, we found the solution t = 1.72 ± 0.30 Gyr, a = 
1.75±0.40, Vi = MiSkms- 1 , Yj = 0.290±0.010 and {Z/X\ = 0.0234±0.0015. 
The position and the evolutionary track of this model (denoted model Mia 
in the following) in the HR diagram are shown in Fig. 1 and 5. The char- 
acteristics of this model are reported in Table 2. Note that the confidence 
limits of each modeling parameter given in Table 2 are estimated as the max- 
imum/minimum values which fit the observational constraints when the other 
calibration parameters are fixed to their medium value. 

Concerning the asteroseismic features of this model, the theoretical variation 
of the large spacings Az/£ for £ = 0,1,2 with frequency was compared to the 
observations (Fig. 2). Table 2 and Fig. 2 show that the mean large spacing of 
the Mia model is in good agreement with the observed value of 55.5 ±0.5 /iHz. 
Moreover, one can see that the Mia model reproduces the observed variation 
of Az/£ with frequency. However, the dispersion of the observed large spacings 
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Table 2 

Models for Procyon A. The upper part of the table gives the observational con- 
straints used for the calibration. The middle part of the table presents the modeling 
parameters with their confidence limits, while the bottom part presents the global 
parameters of the star. 
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around the theoretical curves is slightly greater than expected taking an uncer- 
tainty of 0.57 /zHz, half of the time resolution, on the frequency determination. 
The comparison of the theoretical and observed values of the small spacing 
5uq 2 = v n +i,e=o ~ v n,i=2 between t = and £ = 2 modes is given in Fig. 6. The 
mean small spacing of the Mia model and the theoretical variation of this 
spacing with frequency are compatible with the observed values. However, 
Fig. 6 clearly shows that the observed small spacings are unfortunately not 
accurate enough to provide strong constraints to stellar models. This is espe- 
cially true for the couple of points near 799 and 1136 fxRz that are either split 
by rotation or the secondary peaks are due to noise (see ECBB04). Finally, we 
compare the theoretical p-mode frequencies of the Mia model to the observed 
ones by plotting the echelle diagram (Fig. 8). In this figure, the systematic dif- 
ference (D v ) (see Sect. 3.2.2 of Eggenberger et al. 2004a) between theoretical 
and observed frequencies has been taken into account. Indeed, a linear shift 
of a few /iHz between theoretical and observational frequencies is perfectly 
acceptable, due to the fact that the exact values of the frequencies depend on 
the details of the star's atmosphere where the pulsation is non-adiabatic. 

Model Mia corresponds to the best solution which minimizes the x 2 functional 
of Eq. 1. However, the minimization shows that many models with different 
values of the initial helium abundance and the mixing-length parameter a 
have a x 2 only slightly larger than the one of the Mia model, and constitute 
therefore also good models of Procyon A. This is due to the fact that a de- 
crease/increase of the initial helium abundance Yj can be compensated by an 
increase/decrease of the mixing-length parameter a to match the observed 
position of Procyon A in the HR diagram (see Fig. 3). Thus, we obtain a 
series of models with approximately the same non-asteroseismic features as 
those of the Mia model. Moreover, the mean large spacing of these models 
(which is the only asteroseismic quantity included in our x 2 functional) is 
also very close to the value of the Mia model, since it mainly depends on the 
star's mean density and hence on its radius given that the models have the 
same mass of 1.497 M . This explains why there is a series of models which 
well reproduced the global stellar parameters considered in our minimization. 
Amongst these models, the Mia model is in slightly better agreement with 
these observables. Moreover, when individual theoretical and observational as- 
teroseismic frequencies are compared, one finds that the Mia model is also in 
slightly better agreement with the asteroseismic observations than the other 
models (see Fig. 4). However, given the limited accuracy of the observations, 
and especially of the observed small spacings (which are needed to differen- 
tiate models with the same position in the HR diagram but different ages), 
these differences cannot be considered as really significant. 

To investigate the effects of rotation on the structure of the models and there- 
fore on their p-mode frequencies, we decided to redo the entire calibration 
without including the effects of rotation and atomic diffusion. By perform- 
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Fig. 2. Large spacing versus frequency for the models with a mass of 1.497 M . 
The Mia model includes rotation and atomic diffusion, while the Mlb model has 
been computed without rotation and diffusion (see Table 2). The dots indicate the 
observed values of the large spacing with an uncertainty on individual frequencies 
estimated to half the time resolution. 

ing the same minimization on this grid of standard models, we found the 
solution t = 1.77 ± 0.30 Gyr, a = 1.60 ± 0.35, Yj = 0.280 ± 0.010 and 
(Z/X)[ = 0.0205 ± 0.0015. The position and the evolutionary track of this 
model (denoted model Mlb in the following) in the HR diagram are shown 
in Fig. 1 and 5. The characteristics of this model are reported in Table 2. 
First, we note that the global parameters of this Mlb model are very similar 
to the ones of the Mia model: same effective temperature and approximately 
same luminosity, radius and mean large spacing as those of the Mia model. 
Moreover, the ages of both models are very similar. Of course, the initial 
chemical composition of the models is different, since the Mlb model does 
not take into account the rotationally induced chemical mixing and atomic 
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Fig. 3. Evolutionary tracks in the HR diagram for three models of 1.497 M & with 
different values of the initial helium abundance Y\ and the convection parameter a. 
All models include atomic diffusion and rotation. This shows that a decrease in the 
initial helium abundance Y\ can be compensated by an increase of the convection 
parameter a in order to reach the same location in the HR diagram. 

diffusion. Moreover, the mixing-length parameter of both models is approxi- 
mately equal to the solar calibrated value determined by using the same input 
physics (ck = 1.75 with rotation and atomic diffusion, and a Q = 1.59 without 
rotation and diffusion). 

Concerning the asteroseismic features of the Mlb model, the comparison be- 
tween observed and theoretical variation of the large spacings Ave is shown 
on Fig. 2. This figure and Table 2 show that the mean large spacing and the 
variation of the large spacings with frequency of the Mlb model are in good 
agreement with the observations. Thus, one can see that there are no sig- 
nificant differences between the large spacings of both Ml models; the mean 
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Fig. 4. Differences between calculated and observed frequencies for the models with 
1.497 Mq including rotation and atomic diffusion with different values of the initial 
helium abundance Yj and the convection parameter a. Filled circles correspond 
to the model with Yj = 0.290 and a = 1.75 (model Mia), while open squares 
and crosses designate models with Yj = 0.270, a = 2.0 and Yj = 0.230, a = 2.5, 
respectively. The systematic shifts (D u ) have been taken into account (see text for 
more details). 

large spacing of the Mlb model is only slightly larger than the one of the Mia 
model, which results from the smaller radius of this model compared to the 
one of the Mia model. 

To investigate the effects of rotation on the p-mode frequencies, it is much 
more interesting to consider the small spacings 5vq2 = ^n+i,e=o — ^n,e=2 which 
are mainly sensitive to the stellar core. Indeed, as already pointed out by Di 
Mauro & Christensen-Dalsgaard (2001) and Provost et al. (2002), these spac- 
ings may enable a test for overshoot if sufficiently precise frequencies could be 
determined. To check this, we also computed a model of Procyon A without 
rotation and atomic diffusion, but with an overshoot of the convective core into 
the radiative zone on a distance d ov = a ov min[H p , r core ]. Following Schaller et 
al. (1992), we fixed the value of a m to 0.2, which is the amount of overshoot- 
ing usually chosen for non-rotating stellar models. The evolutionary track of 
this model is shown on Fig. 5. We find that the values of the small spac- 
ings of stellar models including overshooting are slightly smaller than those 
computed for models without overshooting. This can be seen in Fig. 6, where 
the small spacings of two non-rotating models (one with and the other with- 
out overshooting) with the same averaged large spacing are compared. In the 
same way, it is interesting to note that these differences also exist between the 
non-rotating model without overshooting (model Mlb) and the Mia model 
calculated with rotation and atomic diffusion, but without overshooting. In- 
deed, Fig. 6 shows that the non-rotating model (Mlb) exhibits slightly higher 
values of the small spacings than the rotating one (Mia). Thus, we see that 
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p-mode frequencies indicate that the inclusion of rotation mimics an increase 
of the overshoot parameter. This results from the fact that, for sufficiently low 
initial velocities (t> ini < 300 kms -1 ) rotation acts as an overshoot, extending 
the main sequence tracks towards lower effective temperatures (see Fig. 5). 

Even if these theoretical considerations are interesting to point out, it is ev- 
idently not necessary to say that, from an observational point of view, we 
are very far from the precision required to discuss these subtle differences. As 
mentioned above, we can only conclude that the observed frequency differences 
between t = and £ = 2 modes are compatible with theoretical expectations, 
but are not accurate enough to enable us to test the input physics of our stellar 
models. 

4-2 Models with a mass of 1.465 M 

After having calculated stellar models with the mass of 1.497 M determined 
by Girard et al. (2000), we computed a grid of models with a mass of 1.465 M & . 
This mass was deduced using only the WFPC2 measurements with the paral- 
lax of Girard et al. (2000). Accordingly, the value of the observed linear radius 
remains unchanged (2.067 ± 0.028 Rq). All models include shellular rotation 
and atomic diffusion. 

Using this new value for the mass, we found the following solution: t = 1.89 ± 
0.30 Gyr, a = 1.90±0.30, V { = ^iSkms" 1 , Y { = 0.295±0.015 and {ZjX\ = 
0.0229 ± 0.0015. The position of this model (denoted model M2 hereafter) in 
the HR diagram is shown in Fig. 1, while its characteristics are given in Table 2. 

We note that the modeling parameters of this M2 model are very similar 
to the ones of the Mia model. Indeed, the initial chemical composition of 
both models are approximately identical, with a nearly solar initial metal- 
licity which decreases during the star evolution to reach the observed value 
of [Fe/H] s = —0.05. The mixing-length parameter and the initial velocity are 
only slightly larger for the M2 model than for the Mia. Due to the lower mass, 
the age of the M2 model is somewhat larger than the one of the Ml models. 
Note that the large error of Skms" 1 on the initial velocity is due to the large 
error on the observed surface rotational velocity (relative error of about 35 %), 
while the large error on the mixing-length parameter a simply reflects the fact 
that the models are not very sensitive to a change of this parameter. 

Using this lower mass leads to a model whose global parameters still repro- 
duce the observational constraints. However, the agreement between observed 
and theoretical non-asteroseismic parameters is slightly better for the more 
massive Ml models, than for this model. Concerning the asteroseismic fea- 
tures, one can see that the observed large spacings are well reproduced by the 
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Fig. 5. Evolutionary tracks in the HR diagram for three models of 1.497 Mq. The 
continuous line corresponds to the model including rotation and atomic diffusion but 
no overshooting (model Mia). The dotted line denotes the model without rotation, 
atomic diffusion nor overshooting (model Mlb), while the dashed line corresponds to 
a model computed with overshooting (a ov = 0.2), but without rotation and atomic 
diffusion. The dot shows the location of the Mia model. 



M2 model (see Fig. 7). Finally, the echelle diagram of Fig. 8 shows that the 
theoretical frequencies of the M2 model are compatible with the asteroseismic 
observations, even if the agreement is slightly better for the more massive Ml 
models than for this model. However, the limited accuracy of the seismolog- 
ical measurements (especially for the observed small spacings 5u Q2 ) does not 
enable us to formally reject the M2 model. 
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Fig. 6. Small spacing 5^02 versus frequency for the three models with a mass of 
1.497 Mq. The continuous line corresponds to the model including rotation and 
atomic diffusion but no overshooting (model Mia). The dotted line denotes the 
model without rotation, atomic diffusion nor overshooting (model Mlb), while the 
dashed line corresponds to a model computed without rotation or atomic diffusion, 
but with overshooting (a ov = 0.2). The dots indicate the observed values with an 
uncertainty on individual frequencies estimated to half the time resolution. 

4.3 Models with a mass of 1.423 M & 



Finally, we computed a grid of stellar models including shellular rotation and 
atomic diffusion with a mass of 1.423 M Q . This mass was deduced by using 
only the WFPC2 measurements with the Hipparcos parallax (see Sect. 2.1). 
Accordingly, the observed radius of Procyon A is decreased from 2.067 ± 
0.028 R Q to 2.048 ± 0.025 R G . 

By redoing the calibration with this new mass and new observational con- 
straint on the radius, we found the solution t = 2.18±0.30 Gyr, a = 1.80±0.30, 
Vi = 17 ± 8kms-\ Yi = 0.295 ± 0.020 and {ZjX\ = 0.0231 ± 0.0015. The 
position of this model (denoted model M3 in the following) in the HR diagram 
is shown in Fig. 1, while the characteristics of this model are given in Table 2. 
The modeling parameters of this M3 model are also very similar to the ones 
of the other rotating models, apart from a larger age, which results directly 
from the lower mass of the star. 

One can see from Table 2 and Fig. 1 that the agreement between the theoretical 
and observed global parameters of Procyon A is somewhat better for the more 
massive models (models Ml and M2) than for this M3 model. This is especially 
true for the radius of 2.024 R@ which lies in the lower border of the observed 
box (see Fig. 1). This is due to the fact that the radius of this model has to 
be lower than the radii of the other models in order to compensate the mass 
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Fig. 7. Large spacing versus frequency for the models with masses of 1.465 M Q 
(model M2) and 1.423 M & (model M3). Both models include the effects of rotation 
and atomic diffusion. The dots indicate the observed values of the large spacing 
with an uncertainty on individual frequencies estimated to half the time resolution. 



decrease and to keep a value of the mean large spacing close to the observed 
value of 55.5 fxRz. Besides, one can see on Fig. 2 that the theoretical large 
spacings of the M3 model are in good agreement with the observed ones. The 
echelle diagram of Fig. 8 shows that the p-modes frequencies of the M3 model 
are compatible with the asteroseismic observations, even if the agreement is 
better for the more massive Ml models than for this model. 
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Fig. 8. Echelle diagrams for the four models of Procyon A listed in Table 2. Open 
symbols refer to theoretical frequencies, while the filled circles correspond to the 
observed frequencies (see ECBB04). Open circles are used for modes with £ = 0, 
triangles for £ = 1, squares for £ = 2 and pentagons for £ = 3. The models with a 
mass of 1.497 M© are in better accordance with the asteroseismic data than the less 
massive M2 and M3 models. 

5 Conclusion 



By combining all non-asteroseismic observables with the asteroseimic mea- 
surements of ECBB04, we find that the observed mean large spacing of 55.5 ± 
0.5 /iHz favours a mass of 1.497 M Q for Procyon A. Indeed, the lower masses of 
1.465 and 1.423 M lead to models which are in less good agreement with the 
observational constraints. We also determine the following global parameters 
of Procyon A: an age of t = 1.72 ± 0.30 Gyr, an initial helium mass fraction 
Yi = 0.290 ±0.010, a nearly solar initial metallicity {Z/X) { = 0.0234 ± 0.0015 
and a mixing-length parameter a = 1.75 ± 0.40. 
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Note that the lower value of 53.6 ± 0.5 /zHz reported by Martic et al. (2004) 
favours the lower mass of 1.423 M & (Kervella et al. 2004). This can be imme- 
diately understood by recalling that the mean large spacing is proportional to 
the square root of the star's mean density. Since the radius of Procyon A is 
now precisely determined (thanks to interferometric measurements), a lower 
value of the mean large spacing results in a lower value of the mass. 

We also show that the effects of rotation on the inner structure of the star 
may be revealed by asteroseismic observations, if precise measurements of the 
small spacings between £ = 2 and £ = modes can be determined. Indeed, 
there are subtle differences in the frequency dependence of the small spacing 
between rotating and non-rotating models, which results from changes in the 
size and in the outer border of the star's convective core. As expected for small 
initial velocities, these rotational effects are found to mimic the effects due to 
an overshoot of the convective core into the radiative zone. Unfortunately, 
the asteroseismic observations now available for Procyon A are not accurate 
enough to enable us to discuss these subtle differences. 

We conclude that the combination of all non-asteroseismic measurements 
with existing asteroseismic observations puts important constraints on the 
global parameters of Procyon A, but that more accurate asteroseismic data 
are needed to really test the physics of the stellar models. 
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